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Abstract 



Ch ■ A general shell model formalism for the nonmesonic weak decay of the hy- 

•'-H . pernuclei has been developed. It involves a partial wave expansion of the 

X 

. emitted nucleon waves, preserves naturally the antisymmetrization between 

the escaping particles and the residual core, and contains as a particular case 
the weak A-core coupling formalism. The Extreme Particle-Hole Model and 
the Quasiparticle Tamm-Dancoff Approximation are explicitly worked out. It 
is shown that the nuclear structure manifests itself basically through Pauli 
Principle, and a very simple expression is derived for the neutron and proton 
induced decays rates, r„ and Tp, which does not involve the spectroscopic 
factors. We use the standard strangeness-changing weak AA^ — > NN tran- 
sition potential which comprises the exchange of the complete pseudoscalar 
and vector meson octets {7r,r], K, p,uj, K*), taking into account some impor- 
tant parity violating transition operators that are systematically omitted in 
the literature. The interplay between different mesons in the decay of 



is carefully analyzed. With the commonly used parametrization in the One- 
Meson- Exchange Model (OMEM), the calculated rate Tjvm = r„ + Tp is of 
the order of the free A decay rate = r°) and is consistent with ex- 

periments. Yet, the measurements of r„/p = r^/Fp and of Tp are not well 
accounted for by the theory {T^^^^ < 0.42; T^^ > 0.60 T^). It is suggested that, 
unless additional degrees of freedom are incorporated, the OMEM parameters 
should be radically modified. 

PACS numbers:21.80.+a,21.60.-n,13.75.Ev,25.80.Pw 
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I. INTRODUCTION 



The free A hyperon weak decay (with transition rate r° — 2.50 • 10~^ eV) is radically 
modified in the nuclear environment. First, due to Pauli principle the mesonic decay rate 
Tm = r(A — > Nn) is strongly blocked for A > A. Second, new nonmesonic (NM) decay 
channels AN — > NN become open, where there are no pions in the final state. The cor- 
responding transition rates can be stimulated either by protons, Tp = T(Ap — > np), or by 
neutrons, r„ = r(An — > nn). The ultimate result is that in the mass region above A — 12 
the total hypernuclear weak decay rates Fm + T^vm (Fnm = Tn + Tp) are almost constant 
and close to [1]. 

Because of the practical impossibility of having stable A beams, the NM decays in hy- 
pernuclei offer the best opportunity to examine the AS — —1 nonleptonic weak interaction 
between hadrons. Yet, the major motivation for studying these processes stems from the 
inability of the present theories to account for the measurements, in spite of the huge theo- 
retical effort that has been invested in this issue over several decades [2-26] . More precisely, 
the theoretical models reproduce fairly well the experimental values of the total width Fjvm 
(F^£f — r°), but the ratio T^/p = F„/Fp (0.5 < F^^^ < 2) remains a puzzle. 

In the one meson-exchange model (OMEM), which is very often used to describe the 
hypernuclear AN — > A^A^ decay, it is assumed that the process is triggered via the exchange 
of a virtual meson. The obvious candidate is the one-pion-exchange (OPE) mechanism, and 
following the pioneering investigations of Adams [3] ^, several calculations have been done 
in fC, yielding: f£J^) = F^ and Pi^^^^ = 0.1 - 0.2 [10,12,13,21,24]. 

The importance of the p meson in the weak decay mechanism was first discussed by 
McKellar and Gibson [4]. They found that, because of the sensitivity of the results to the 
unknown AA^p vertex, the estimates for Fatm could vary by a factor of 2 or 3 when the 



^McKellar and Gibson [4] have pointed out that this publication contains a very important error 
and that the decay rates given by Adams [3] should be multipUed by 6.81. 



potential Vp was included. (See also Ref. [5].) The present-day consensus is, however, that 
the effect of the p-meson on both Tnm and Tn/p is small [10,12,13,15]. 

Until recently there have been quite dissimilar opinions regarding the full OMEM, which 
encompasses all pseudoscalar mesons (n^rj^K) and all vector mesons {p,u;,K*). In fact, 
while Dubach et al. [12] claimed that the inclusion of additional exchanges in the tt + p 
model plays a major role in increasing the n/p ratio, Parreno, Ramos and Bennhold [13], 
and Sasaki, Inoue and Oka [21] argued that the overall effect of the heavier mesons on this 
observable was very small. However, the two latter groups have recently corrected their 
calculations for a mistake in including the K and K* mesons, and so their estimates of T^/p 
have augmented quite substantially [22,24]. Almost simultaneously, Oset, Jido and Palomar 
[25] have also shown that the K meson contribution was essential to increase Tn/p. However, 
the experimental data have not been fully explained yet. 

In the last few years, many other attempts have not been particularly successful either 
in accounting for the measured T^/p ratio. To mention just a few of them: 1) analysis 
of the two-nucleon stimulated process KNN NNN [7,9,14], 2) inclusion of interaction 
terms that violate the isospin AT = 1/2 rule [16,20], 3) description of the short range 
baryon-baryon interaction in terms of quark degrees of freedom [17,21], and 4) introduction 
of correlated two-pion exchange potentials besides the OPE [18]. Consistent (though not 
sufficient) increases of the nip ratio were found in the last two works. (For instance, T^/p 
was boosted up to 0.36 for the decay of [18].) In fact, only Jun [26] was able so far to 
reproduce well the Fjvm, Fp and T^/p data. He has employed, in addition to the OPE, an 
entirely phenomenological 4-baryon point interaction for short range interaction, including 
the AT = 3/2 contribution as well, and has conveniently fixed the different model coupling 
constants. Let us also note that after the present work had been completed, Itonaga et 
al. [27] have updated their studies and have performed extensive calculations of the NM 
decays in the mass region A < A < 209, which have revealed that the correlated-27r and lu 
exchange potentials significantly improve the T^/p ratios over the OPE results. 

In the OMEM's, a weak baryon-baryon-meson (BBM) coupling is always combined with 



a strong BBM coupling. The strong one is determined experimentally with some help from 
the SU(3) symmetry, and the involving uncertainties have been copiously discussed in the 
literature [28-31]. It is the weak BBM couplings which could become the largest source 
of errors. In fact, only the weak A^'Att amplitude can be taken from the experiment, at 
the expense of neglecting the off-mass-shell corrections. All other weak BBM couplings 
are derived theoretically by using SU (3) and SU (6)^ symmetries, octet dominance, current 
algebra, PCAC, pole dominance, etc. [6,12,13,32-37]. Assortments of such methods have 
been developed and employed for a long time in weak interaction physics to explain the hy- 
peron nonleptonic decays. Specifically, to obtain the weak BBM couplings for vector mesons 
the SU {6)yj symmetry is used, which is not so well estabhshed as the SU{3) symmetry is. 
Moreover, the results derived by way of the SU (6)^ symmetry depend on the contributions 
of factorizable terms ay and qt, which were only very roughly estimated [12,32,34]. Well 
aware of all these hmitations, McKellar and Gibson [4] have allowed for an arbitrary phase 
between the p and the tt amplitudes in the tt + p model. The same criterion was adopted 
by Takeuchi, Takaki and Bando [5]. 

We wish to restate that the OMEM transition potential is purely phenomenological 
and that it is not derived from a fundamental underlying form, as happens for instance in 
the case of electro-magnetic transitions or the semileptonic weak decays. Only the OPE 
model is a natural and simple extrapolation of the mesonic decay mechanism of the A to 
the NM process: the weak BBM coupling is identical to that used in the phenomenological 
description of the free A, and the strong vertex is the one traditionally used in describing the 
ttNN vertex. The assumption is that this is a valid approximation, although the pion is off 
the mass shell. Accordingly, all modern interpretations of the NM weak decay use the OPE 
as the basic building block for the medium and long range part of the decay interaction. On 
the contrary, the full OMEM is not used very often and, in place of the one-meson 77, • • • 
exchanges, other mechanisms are employed as refer to above. One should also keep in mind 
that both the strong and weak BBM couplings, as well the meson masses, can become 
significantly renormalized by the nuclear environment [38]. 



The high momentum transfer in the NM decay makes the corresponding transition am- 
phtude very sensitive to the short range behavior of the NN' and NA interactions. In fact, 
quite recently it has been pointed out that the final state interactions (FSI) have very large 
influence on the total and partial decay rates [24] (see also Ref. [27]). ^ Due to the same 
reason one could expect that the nuclear structure effects not included in the main field 
(such as the RPA or pairing correlations, higher order seniority excitations in the initial and 
final states, etc) should not play an important role. Yet, it could be useful to understand 
this issue more genuinely and to get a more complete control on the nuclear structure aspect 
of the problem. These are the main motivations for the present work. 

The only existing shell model framework for the hypernuclear decay is the one based on 
the weak-coupling model (WCM) between the hyperon and the (^4 — 1) core [8,13,18,27]. It 
involves the technique of coefficients of fractional parentage, and the spectroscopic factors 
(SF) explicitly appear in the expressions for the transition rates. Yet, in nuclear structure 
calculations it is in general simpler to evaluate the transition probabilities directly from 
the wave functions, instead of doing it via the SF. Here we first develop a fully general 
shell model formalism and then we work out thoroughly the extreme particle-hole model 
(EPHM) and the Quasiparticle Tamm-Dancoff Approximation (QTDA) for the even-mass 
hypernuclei. 

Owing to the above mentioned characteristics of the OMEM it might be legitimate to 
ask whether it is possible to account for all three data Ftvm, Fp and F„/p by not fully 
complying with the constraints imposed by the SU{3), SU{6)ui and chiral symmetries on 
the BBM couplings. To find out in which way should these parameters be varied we perform 



^The FSI also make hard the extraction of the n/p ratio from the experimental data [39-43], and 
to surmount this difficulty Hashimoto et al. [43] have quite recently combined the Monte Carlo FSI 
internuclear cascade models from Ref. [14] with the geometry of the detectors. Moreover, Golak et 
al. [47] have shown that the FSI, in principle, hinder the measurement of the n/p ratio in ^H. 



a multipole expansion of the transition rate in the framework of the EPHM, which unravels 
in an analytic way the interplay between different mesons in each multipole channel. 

Attention will be given also to the parity violating potential, since there are several 
typographical errors in the recent papers [13,19,24], regarding this part of the transition 
potential. We will also consider some important contributions due to the vector mesons, 
which, although always included in the description of the nuclear parity violation [32,44-46], 
have been so far neglected in all studies of the NM hypernuclear decays, except those of 
Dubach et al. [12,34]. 

The outline of this paper is as follows: The general shell model formalism for the hyper- 
nuclear KN — > NN weak decay is developed in Sec. II. The nonrelativistic approximation 
for the effective Hamiltonian is presented in Sec. III. The EPHM and QTDA are explained 
in Sec. IV, where the multipole expansion of Fjvm is also done. Numerical evaluations of 
the ^^C — > ^^C + nn, and ^^C — > ^^B + pn decay rates, are carried out in Sec. V, and 
the conclusions are presented in Sec. VI. The formulae for the nuclear matrix elements are 
summarized in the Appendix. 



II. TRANSITION RATE 



The decay rate, of a hypernucleus (with spin Jj and energy £j) to residual nuclei (with 
spins Jp and energies Ep) and two free nucleons (with total spin S and energies Cp and ep), 
follows from Fermi's Golden Rule 



r = 27r f\{pPSMs;JFMF\V\JiMi)\''S{ep + ep + SF-£i)j^^. (2.1) 

Here, V is the weak hypernuclear potential, the wave functions for the kets \pPSMs; JpMp) 
and I J7M7) are assumed to be antisymmetrized and normalized, and a transformation to 
the relative and center of mass (cm.) momenta, p and P, is already implied, i.e., 

P = ^(Pl-P2), P = Pl+P2. (2.2) 

It is convenient to define the quantity 

I{p,P) = (An)-^ J2 j dVtpdnp\{pPSMs]JFMF\V\JiMi)\^, (2.3) 

SMsJpMp ■' 

and rewrite (2.1) as: 

r = ^ / de^e{Ap - 6)J(p, P), (2.4) 

where P — 2^/M^, p — Mf^{A — e), Ap — Si — Sp — 2M]y, and Mjy is the nucleon mass. 

The partial wave expansion of the wave function of the non-antisymmetrized two-particle 
ket \PpSMs) is then performed: 

{rRs,S2\PpSMs) ^ {4nf ^ z'+''Yi:^{p)Y*^{P){rKs,S2\plm, PLM, SMs), (2.5) 

ImLM 

where 

{rRsiS2\plm, PLM, SMg) = Yim{r)YLM{B:)ji{pr)jL{PR)xsMs{siS2), (2.6) 
describe the spherical free waves for the outgoing particles, 

r = ri-r2, R=^(ri + r2). (2.7) 



are the relative and cm. coordinates, and I and L are the quantum numbers for the relative 
(1) and cm. (L) orbital angular momenta. After performing the angular integration in (2.3) 
we obtain: 

I{p,P)= Yl E \{plm,PLM,SMs;JFMF\V\JiMi)\^, (2.8) 

SMsJfMf ImLM 

which goes into 

I{p,P)^ \{pPILXSJ,Jf;JiMi\V\JiMi)\'^. (2.9) 

SlLXJJp 

when the angular momentum couplings: 1+L = A, A+S = J are carried out. The quantum 
number Mj is superfluous and will be omitted from now on. 
The transition potential is written in the Fock space as: 

Y {pPlLXSJ\V\uj^J) • (2-10) 

where, in the same way as in (2.1), a transformation to the relative and cm. momenta is 
implied. Here, ja and are the single-particle shell-model states of the decaying particles, 
and = (-)^+"^a,_^ [49]. One gets 



SILXJJ^ 



Y {pPlLXSJ\V\j^jr,J){Ji\\ {alal)j\m 



JNJA 



(2.11) 



where the transition amplitudes (J/|| (oj^v'^L) j 11"^^) reduced with respect to the angular 
momenta, the label a stands for different final states with the same spin Jp, and J = 
V2J+1. 

The effective weak hypernuclear interaction is isospin dependent, i.e., 



V{r,SASN:tAtN)^ Y^r{r,SASN)Tr, % ^ I , (2.12) 

and therefore the nuclear matrix elements have to be evaluated in the isospin formalism. 
This implies that (2.11) goes into 



SXILTJJ'^ 



JNJA 



(2.13) 



where 



X 



Y,{pPlLXSJ\Vr\jAjNJ){TMT = mt^ + mtMrnt^rritr,), (2-14) 



is the antisymmetrized nuclear matrix element, and rritj, — \ and rrit^ = rut^ 



-\. It 



is assumed, as usual [13], that A behaves as a "2) isospin state. In that way the phe- 
nomenological AT — | rule is incorporated in the effective interaction. Note that in (2.13) 
and (2.14) mt^ = Mt - rrit^. 

To evaluate {pPlLXSJ\VT-\jAjNJ) one has to carry out the jj — LS recoupling and the 
Moshinsky transformation [50] on the ket IJaJnJ) to get 



(pPlLXSJ\Vr\jAjNJ)=jAjN E ^'S' 

A'5'nlNL 



5 JA 



Ir ' 



2 

N \ Jn 

X' S' J 

X {nmiX'\nAlAnNlNy){pPlLXSJ\Vr\nmiX'S'J), 



(2.15) 



where (•••!•••) are the Moshinsky brackets [50]. Here, I and L stand for the quantum 
numbers of the relative and cm. orbital angular momenta in the AN system. The explicit 
expressions for the transition potentials are given in the next section, and the formulae that 
are needed to evaluate the matrix elements (pPiLA5'J|Vr|nlNLA'S" J) and {T MT\%-\mtj^mtj^) 
are summarized in the Appendix. 

When the hyperon is assumed to be weakly coupled to the A — 1 core, which implies 
that the interaction of A with core nucleons is disregarded, one has that | J/) = \{JcjA)Ji), 
where Jc is the spin of the core. Prom 



Jc Jl JA 

J jjv Jf 



{JcWaLmJW, (2.16) 



we obtain 



SlLXJTJf 



X 



J2M{pPlLXSJT;jAjN:mt^) < 



' {Jc\\a]^mJ\JF) 



(2.17) 



Occasionally it could be convenient to include the isospin coupling as well into 



{JcTcMTa\W.^J\JFTFMT, 



(2.18) 



[Jn\\a\.,^. \\Jf)^ and work with the spin-isospin reduced parentage coefficients 

jt I 

{JcTc\\\a^^.J\\JFTF) - Tc ^TFMrAmjTcMr,) 

where Tc, Mt^ and Tp, are the isospin quantum numbers of the core and residual nuclei, 
respectively. In this case 



JfTpSlLXJT ^ 



X 



E M {pPlLXS J ; jAjN , mtj, ) 



JN 



Jc Ji 3k 

J Jn Jf 



{JcTc\\\al,J\\J^TF) 



(2.19) 



Thus, knowing the transition potential V and the initial and final nuclear wave func- 
tions I J/) and \Jf) (or \Jc) and | Jf)), we can evaluate the transition rate (2.4), with the 
integrations going up to 



Ma-Mn + + e 



3K ^ ^ONrntj^ 1 



(2.20) 



where e^^ and ej^^^^ are the single particle energies. 



III. EFFECTIVE INTERACTION 



As the reduction of the relativistic one-meson exchange t-matrix, to the nonrelativistic 
effective potential V, is in the hterature [4,6,10-13,24,34,37] it will not be repeated here. 
For the parity conserving (PC) potential we will just list a few results that are indispensable 
for establishing the notation and for the final discussion. More attention will be given 
to the parity violating (PV) potentials. In dealing with them some tricky details appear 
concerning the passage from the momentum space to the coordinate space. We first illustrate 
the procedure for one pseudoscalar meson (tt) and one vector meson (p) , and afterwards we 
generalize the results to all six mesons. 

The effective strong (S) and weak (W) Hamiltonians read 



(3.1) 



where Gpm^ is the weak coupling constant, and are the baryon fields, tt and p are 
the meson fields, r is the isospin operator, M the nuclcon mass, and M the average between 



the nucleon and A masses. The isospin spurion 



is included in order to enforce the 



empirical AT = \ rule [13]. 

The corresponding nonrelativistic t-matrix in the momentum space (with the hyperon A 
being always in the first vertex) is: 



^7r(q) = -Ta ■ Tn- 



At,{(Tn ■ q) + B^(o-A • q)(<TAr • q) 



ml + q2 



^, Up((TA X (Tat) •q-2^;,<TA-Q + Bp(<TA X q)((TAr X q) -B;, 

q, Q = -'^A • 'Tjv — ^- — ^ ^, 3.2 

where the couphng constants Am-, A'j^, Bm and are defined in Table I, and 

q = p'-p; Q = ^(p+p), (3.3) 

with p' and p being, respectively, the relative momenta for the initial and final states. (We 
have adopted this labeling to be consistent with (2.2).) In the momentum space the potential 
reads: 

(piP2|l^|p'iP2) = -(27r)^5(p; + p'2 - pi - P2)t(q, Q), (3.4) 
and in order to arrive to the coordinate space the Fourier transform is applied: 



X exp{i[p'i • r'l + P2 • r2 - pi • ri - p2 • r2]}. (3.5) 
After some trivial integrations and the coordinate transformation: 

x = r-r'; X = i(r' + r), (3.6) 

we get 

(rir2|y|r;r^) = -5(R'-R) / exp[i(Q • x + q • X)]t(q, Q). (3.7) 

To carry out the integration on q and Q we make use of the result: 



/ 



/ 
/ 



[Q ,,.2 , -9 = o^(r - r)fM (n, 



(27r)3(27r)3^ m|^ + q2 2 
r/O (Vn „i(Qx+qX) 
^ {a,- q)(a2 • q) , = -6{r' - r)[/^(r)(cT, ■ a^) + /^(r)5i2(f)], 



(27r)3(27r)3 



m|f + q2 



ciQ dq , i(Q.X+q.X) 



-5(r'-r)[2/j^(r)(cTi-cr2)-/^(r)5i2(f)], 



where 



(3.8) 



'S'i2(r) = 3(<Ti • r)(o-2 • f) - (cti • 0-2) 



124:71 



Y2{f) ■ [<Ti X <T2]2, 



(3.9) 



is the tensor operator, and the radial dependence is contained in: 



fl7^(r) = [V, /M(r)] = V/M(r) = r^/M(r) = r/;,(r), 



fl^^ W = {V, /M(r)} = V/M(r) + 2/M(r) V, 



/mW = o m^/M(r)-(5(r) 



m 



M 



3 3 
1 + + 



(3.10) 



with V = V12 = Vi = — V2, and 



/M(r) 



-; r=|ri-r2|, 



Thus (3.7) reads 



V mMfJ 



(3.11) 



(rir2|y|r;r'2) = 6{r' - r)6{R' - R)V{r), 



(3.12) 



where the transition potential for the n + p model is: 



+ BMr) - iA^a^ ■ f(-)(r) - • f(+)(r) + A,{<t^ x a^) ■ f(-)(r)} 



(3.13) 



The complete potential can now be cast in the form (2.12), with the isoscalar (r;, a;) 
and isovector (tt, p) mesons giving rise to Vo and Vi, respectively, while the strange mesons 
{K,K*) contribute to both. We get: 



+ icrA 



+ (o-A X (Tn) ■ [AJt\r) + Ax^ij 



K* \ 



+ (o-A X o-at) ■ \Apflr\r)+AKii^KKr) 



(3.14) 



for the PV potential, and 



+ SAN{r) B,f^{r) + BkJUt) - B^f^{r) - Bk^J^t) 



+ BlUr)+B'fK^ir), 



+ ^Aiv(r) [B./J(r) + BkJUt) - B,fJ{r) - BK^fUr) 



+ B'^f,ir)+B'^.fK*ir), 



(3.15) 



for the PC potential. The overall coupling constants Am, A'm, and , B^, are listed in 
Table I, with the weak couplings for kaons defined as: 



\-Uk ^ 



2 ' 



R - _L nPC. R _ 



iPcy 



PC,V 



2 



2 ' 



5 



C 



PC,T 



2 



c 



PC,T 



K* 



(3.16) 



The C's and D's are given in Ref. [13]. The operators that have been habitually omitted in 
V^^(r, S1S2) are those that are proportional to A!^- 



IV. NUCLEAR MODELS AND MULTIPOLE EXPANSION 
A. Extreme Particle-Hole Model 




J, J 

FIG. 1. Diagramatic representation of the hypernuclear NM weak decay, from the Iplh state 
\jAja^',Ji) to the 2h state Ija^j^^Jp), while two nucleons with momenta pi and p2 are emitted 
into the continuum. S and W are the strong and the weak vertices, respectively, and M is a 
non-strange meson. 

The simplest nuclear shell model is the EPHM, in which: the hypernucleus fZ is de- 
scribed as a A-hyperon in the single particle state |ja) and a hole state \ja~^) relative to 
the core, while the residual ^~^Z and — 1) nuclei are represented by the two hole 

states \ja~^jb~^) with respect to the same core. As illustrated in Fig. 1 | J/) — > IjaJ^''^; Ji), 
\Jf) \ja~^jb~^', Jf) and \ Jc) |ja^^)- The parentage coefficients in (2.17) read 

(JCIISUJIJ^) = (_)^-+^a+^yi + 5„,J^. (4.1) 

In particular, for the initial state is 

J, = -^^^ = \^sy2K lpz/2n~\ 1), (4.2) 

and the final states are: 



|(lp3/2n"i)2;0,2),|lp3/2n-Hsi/2n-i;l,2), 



(4.3) 



|(lp3/2n Hp3/2P ^);0, 1,2,3), |lp3/2n Hsi/ap ^1,2), 

for An — > nn and Ap — > np, respectively. Here |0) is the ^^C particle vacuum. As there is 
only one hole state for each parity, the parentage coefficients with different — do not 
interfere among themselves. After summing up on the final states the integrand (2.17) can 
be cast in the form 



ImAP,n = T.KAPy2) E M\pPlLXSJT;j^j,,mt,), 



(4.4) 



SXILT 



where 



ja+jb 



Jmt^mta 



dp 



JF = \3a-3b\ 



ja jb Jf 

J Ji Ja 



(4.5) 



are geometrical factors which come from the Pauli principle. Their explicit values for the 
lsi/2, lp3/2, lpi/2 are listed in Table II. 



B. Beyond Extreme Particle-Hole Model 

The EPHM can be straightforwardly improved by going to the quasiparticle represen- 
tation. In fact, for all even-mass hypernuclei the initial and final states can be expressed 
as: 

3A3a 

\JF) = T.ChbjA^i^i)j,\^^^)^ (4.6) 

3a3b 

where 6^ = UjOj — Vjaj is the quasiparticle creation operator [48], \BCS) is the BCS vacuum, 
ja is always a neutron state, while jb can be both a neutron and a proton orbital. Note that 
because of the lack of hyperon-hole states, the backward going RPA contributions do not 
appear and one has to work within the QTDA. Prom (2.13) we get 



ILXSJJ^T 



I I ia Jl 7A 

3A3a3b I J jl, Jp 



2 



(4.7) 



The residual interaction in the final nuclei redistributes the transition rates among the states 
with the same spin and parity. But, as the NM decay is an inclusive process, i.e., the partial 
transition rates are summed up coherently over all final states, such a rearrangement plays 
only a very minor role on the total rates. (The same happens, for instance, in the neutrino- 
nucleus reactions and in the meson capture [51].) Therefore, it is justifiable to approximate 
the final wave functions by their unperturbed forms, i.e., \Jp) = {b]h]^j \BCS) and 
^fajbJp ^ ^o'Jajb- addition, one assumes that the hyperon is always in the lowest lsi/2 

state the last equation takes the form of (4.4), i.e., 

^m^MP) = E K,,jiJa)Cl,y,, E M\pPlLXSJT;uj,,m,,), (4.8) 

jaJbJ SXILT 

Only the orbitals lsi/2, 1^3/2 and lpi/2 will be used. In this case, as seen from Table II, 
PpjiPi/2) — -f'pj(P3/2), which implies that in the case of protons the summation on ja can be 
performed analytically. Thus, as J2ja ^jajA ~ finds out that Tp does not depend at all 

on the initial wave function. From the same table one also finds out that -F^j(ja) = FpjUa), 
except when — ja. So, one can expect as well only a weak dependence of r„ on | J/). This 
fact is verified numerically later on. 

In summary, we end up with a very simple result for the transition rates: 

Tm*, = E^i^^;J0a)7^^4J, (4.9) 

3b J 



where 

16MI f^.b 

SILXT 

Clearly, the EPHM is contained in (4.9) with the occupation numbers Vj,^ equal to one for 
the occupied states and to zero for the empty states. 



^t,. = / de^e{A,^ - e) M\pPlLXSJT; j^,,m,,). (4.10) 



C. Multipole Expansion 



The EPHM is particularly suitable for performing the multipole expansion of the in- 
tegrands Im^- Thus we carry out both the Racah algebra in (2.15) and the summations 
indicated in (4.4), keeping in mind that the allowed quantum numbers {IL} are: {00} for the 
Si/2 state, and {01} and {10} for the P3/2 state. To simphfy the results we take advantage 
of the relations 

1/4 



(P0|10) = (0 b'/'e-^^'^'/\ 



(P1|11) = -L(0 6V2pe-(^6)V4^ (4.11) 



for the radial integrals (PL|NL) defined in (A2), and introduce the ratio 



(Pl|ll) 



2 



(4.12) 



(P0|10)_ 

which allows us to work only with the L = overlap (P0|10). Thus, from now on the label 
L will be disregarded, and to identify the Si/2 and P3/2 pieces of the I = strength we will 
use the ratio R, which appears only in the last term of (4.4). The results of the multipole 
expansion for both PC and PV potentials are displayed below. 

1. Parity conserving contributions 

The matrix elements of the PC operators fM{r), f^{r){cr\-cr]sf), and /M(^)'S'AAr(?^), given 
by (Al), can be expressed by means of the radial matrix elements (A2) and (A3), or more 
precisely through the moments 

c^(p,p) = s;,(pI|/m|ii)(po|io), 

1 for TT, T], K 



S^(p,P) = BM(pl|/j^|ll)(P0|10)x 



2 for p,uj,K* 



T^(p,P) = BM(pi|/^|ll)(P0|10)x <( ^ ''''''^ . (4.13) 

—1 for p,u},K* 



Introducing the notation: 







Co = Co, + Cko 
So = + S^^ + Sko + Si^* 
To = + T^^ + + Tx* 



Ci — Cp + Ck* , 

Si = + Sp + Ski + Sk*, 



for the isoscalar (r = 0) and the isovector (r = 1) matrix elements, one gets: 



2p = 2 (1 + i?) 



3 (S°)' + 9 (S?)' + (C°)' + 7 (C?)' + 6 (3Tf - T^^^) 



- 4C°C° + 12C°S? - 6C°S? - 6C$S° 
+ 6 [Slf + 42 (S})' - 24SjSi + 2 (cj)' + 6 (C})' - 24CiSi 



for the decay Ap — > np, and 

7R 



In=[l + f) (3SU3S?-C°-C?)%H(si + q + S} + C}) 



for the decay An — > nn. 



2. Parity violating contributions 
The PV matrix elements (A5) are reduced to the nuclear moments 

p'i,(p,p)=^M(p/i/yiii)(poiio), 

Qi{,(p,P) = ^'^(pi|/l+)|ll)(P0|10), 
where the radial integrals are defined in (A8). Using the notation, 

P?? = P»j — Qk*, ^Ko = P/fo ~ ^u), PftT* = Pk* + Pa;, 
Ptt = Ptt — QjiT*, Pjfi = — Qp, Pr* — Pr; + Pp, 



we obtain: 



X, = 2 (1 + i?) [3 {pi'Y + {pfY + 3 {p^y + {PiY + 10 {p'^,y + 2 {p^^,y 

OplOplO I oplO /'oplO plO I /1 510 OPIO \ I /iplO /QplO pio \ /iplO plO 

+ 14 (P-)^ + 2 (Pjf + 8 (P- )^ + I (P- + 14 (P- )^ + ^-^ (P- " 



3 V "0/ ■ V "1/ ■ 3 
+ 4P?P^\ - 4Pf (2P^. + 2P^\ - P^„ + 4P^^ - 2P^.) + 4Pi\ {Pf - P^^ + P^,) 

I /|p21 /'9p21 p21 \ I 4~2i ~2i 

+ I (P^o)' + 6 (P- )^ + I (P- )^ + 6 (P- )^ - 4P- P- 



"^pOl (pOl qpOl I QpOl \ I /ipOl /pOl QpOl \ 



(4.19) 



for the Ap — > np decay, and 

43i?\ 

1 
R 

9" 

14i?\ 



pio , plO\ /plO , plO ' 



j K + p^°0 (p^v p;° + p^„ + p],° ) 

+ ^(p^^ + pf + p^o + p^\)' 

+ 2 (p- + p- ) (p- + p- + pf + pf + p- + p 
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(4.20) 



for the An — > decay. 



V. NUMERICAL RESULTS AND DISCUSSION 



The numerical values of the parameters, defined in Table I and necessary to specify 
the transition potential, are summarized in Table III. For the sake of comparison all cutoffs 
appearing in (5.1), as well as all coupling constants, were taken from Ref. [13], where, in turn, 
the strong couphngs have been taken from Refs. [28,29] and the weak ones from Ref. [12]. 
The energy difference Aj^mt^ in (2.20) is evaluated from the experimental single nucleon 
and hyperon energies, quoted in Ref. [8]. 



The finite nucleon size (FNS) effects at the interaction vertices are gauged by the 
)nopole form factor FjJ^^''(q^^ 
and (3.11) must be replaced by 



monopole form factor FjJ^^''(q^) = which implies that the propagators in (3.10) 



/M(r) ^ /M(r) = /M(r) - /A^(r) - ^^^^%-^/A^(r), 

fnir) - lUr) = /^(r) - fijr) - ^(A^, - ml){AMr - 2)U^{r), 

fLir) - flir) = fLir) - fijr) - ^(A^ - mL)(AMr + 1)/am 

fuir) - 7M(r) = fuir) - A^(r) + ^^^^^^^/Am W, (5.1) 



where fhf^^ij-) has the same structure as fuij') but with uim — ^ Am- 

The initial and final short range correlations (SRC) are taken into account, respectively, 
via the correlation functions [13], 

9f{r) = l-Jo{qcr), (5.2) 

with a = 0.5 fm, /? = 0.25 fm'^ and 7 = 1.28 fm, and = 3.93 fm-^ 

It is a general belief nowadays that, in any realistic evaluation of the hypernuclear NM 
decay, the FNS and SRC have to be included simultaneously. Therefore, in the present paper 
we will discuss only the numerical results, in which both of these renormalization effects are 
considered. Under these circumstances, and because of the relative smallness of pion mass, 
the transition is dominated by the OPE [13]. 



The major part of the numerical calculations were done in the EPHM where the only free 
parameter is the harmonic oscillator length b. The most commonly used estimate is 6 = A^^^ 
fm [49,48], which corresponds to the oscillator energy gives b — 1.51 

fm. For light nuclei it is sometimes preferred to employ hu — 45 — 25A~^/^ MeV, which 
yields b = 1.70 fm. Moreover, a y\ particle in a hypernucleus is typically less bound than 
the corresponding nucleon and hence 6a could be larger than b^. For instance, in Ref. [13] 
was used 6 = (6a + 6Ar)/2 = 1.75 fm, which comes from bN — 1-64 fm and 6a = 1-87 fm. 
As there is no deep motivation for preferring one particular value of b, the numerical results 
will be exhibited for both 6 = 1.51 fm and b = 1.75 fm. 

First, a few illustrative results, obtained in the EPHM (Eq. (4.4)) and the simphfied 
version of the QTDA (Eq. (4.8)), are displayed in the Table IV. The hyperon- nucleon 
interaction in the later approach was taken to be a simple S force, which has been recently 
used with success as the nucleon-nucleon interaction to explain the weak decay processes 
in ^^C [51]. The resulting pairing BCS factors were: Vg^^^ — 0.9868, Vp^^^ = 0.8978 and 
~ 0.6439. Although we have expected to obtain small differences between the EPHM 
and the QTDA, it came as a surprise that they turned out to be so tiny. Thus, henceforth 
only the first one will be used. 

Next, we combine results from Table IV with the multipole expansion done in the previous 
section to find out the roles played by different mesons. Note that the formulas (4.19) and 
(4.20) depend on the ratio (4.12), and it was found numerically that the approximation 

i?=l, (5.3) 

reproduces fairly well the exact calculations. This estimate helped us to formulate the 
following comments: 

PC potential; The dominant contributions to Fp and F„ come from the I = matrix 
elements, while the I = 1 wave contributes relatively little: = 2% to Fp and = 10% to F„. 
On the other hand, for the parametrization displayed in Table III, one finds that: (1) the lu 
and K* mesons mainly cancel out in Cq, as the p and K* mesons do in C?, and (2) the matrix 



elements 5° and are small in comparison with (ST^'^ — Tg*^), which makes Tp to be large 
vis-a-vis r„. Thus, using the estimate (5.3), one ends up with the following approximate 
result for the PC contributions: 

Ip + I^^Ip^ 2i;'/' ^ 24 (3Tf - T^o)' . (5.4) 

From Table III and (4.14) one can also see that: (i) the lu and K* mesons contribute 
coherently with the pion, while the remaining three mesons contribute out of phase, and 
(ii) the different vector meson contributions have the tendency to cancel among themselves. 
As is shown in Table IV, the overall effect is a reduction of the pion transition rate by 
approximately a factor of two. 

PV potential; As in the PC case, the dominant PV transition strengths come from 
the I = wave, through the P]^ moments. The I = 1 wave from the P3/2 state gives rise to 
I — and I — 2 outgoing channels. The first one can always be neglected, while the second 
one contributes with 15% to Fp and with 2% to F„, when only the tt meson is considered. 
After including all mesons these percentages drop to 6% and 1%, respectively. Also here 
the partial Si/2 and P2/2 contributions are approximately equal for all mesons in the proton 
induced channel and notably different in the neutron induced channel. 

From Table IV it is easily found that the most important PV contributions arise from the 
moment, and from its interference with the P^°, P)^ and Pj^. moments. Thus, retaining 
only the most relevant terms in (4.19) and (4.20), the following rough estimates are obtained: 



Ip ^ 2Ip'^^ ^ 4 



3(Plf + 3(P^°)%(P],«„)^ 
+8Pf {2P'l - PZ + 4P1°. - 2P]^, - 3Q]°.) , (5.5) 



and 



97 

"~ 18 



■'■^pio fplO , plO 
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78 

These relations are notably more comphcated than (5.4). Nevertheless, it can be concluded 



that: (1) the t] meson is only significant for r„, and (2) the K and K* mesons increase both 
transition rates, but in a different way. 
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FIG. 2. Matrix elements of the radial operators /^., /^*, fK*^ ^^'^ Ik*^ ^ function of the 
energy. 

Before proceeding it is worth to say a few words on the "new" nuclear moments Q^j and 
compare them with the well known moments Pj^j. As seen from (4.17) and (A8)-(A10) they 
basically differ in the radial dependence. Specifically, we discuss the radial matrix element 

(p, l|/i;J|10) = {p, llfkAW) + {p, 1|/^.|10), (5.7) 

which appears in Q^*, together with the usual matrix element 

(p, 11/^110) = (p,l 1/^*1 10) - (p,l|/|,|10), (5.8) 

which is contained in Pj^. . The overline indicates that both the FNS and SRC are included, 
as explained in the Appendix. 

As can be seen from Fig. 2, the matrix elements of /^* and have opposite signs, 
and as a consequence the matrix element of f^J is larger in magnitude than that of f^j^*- 
A rough approximation for the mean values is: 



As Ak* ^ A'j^* (see Table III) we end up with the estimate: 

(5.10) 

Thus, equations (5.5) and (5.6) show that the K* meson mainly contributes through the 
moments P^*, augmenting the magnitude of F^^ and diminishing that of F^^. The matrix 
elements Q^*, in contrast, reduce both transition rates. 

Furthermore, the equation (4.18) indicates that each vector moment Q^^^^ k* accom- 
panied by a pseudo-scalar moment Pi-°r),if- Both integrals (pl|/^''|10) are negative for all 
mesons. Then, using the values of the coupling constants An,n,K and A'p ,^ x* listed in Table 
III, it can be inferred that Qp^L.i^* ^.nd P^^^jj^k moments mostly add incoherently. 

The experimental results for the total transition rate T^m, the proton partial width 
Fp, and the ratio F„/p in are displayed in Table V. In the same table the theoretical 
estimates are also shown, grouped as: 

1. Calculation A. All the parametrization is taken from Table III, and the following cases 
are shown and commented: 

• (tt): The simple OPE model accounts for T^m, but it badly fails regarding Fp 
and Tn/p. 

• (PS): When 77 and K mesons are included, the total transition rate is only slightly 
modified, while Fp and F„/p change significantly, coming somewhat closer to the 
measured values. 

• {PS + K*): The incorporation of the K* meson increases Fjvm and Fp, decreases 
F„/p, and in this way worsens the agreement with the data. 

• {PS + V) : The results are not drastically modified when all vector mesons are 
built-in. 



Q 



0.3 



pio 



• {PS + V{P)): All 6 mesons are included, but only the PV moments Pm are con- 
sidered. The importance of the new moments Qm is evident from the comparison 
with the previous case. 

The main conclusion is that it is not possible to reproduce simultaneously the data for 
all three observables Fjvm, Tp and r„/p, when the BBM couphng are constrained by 
the SU (3) and SU (6)^ symmetries. 

2. Calculation B. We discuss now what happens when the just mentioned constraints are 
relaxed, and the FNS and SRC parametrizations, as well as the the pion couplings, 
are kept unchanging. That is, the transition potential is considered to be given by 
a series of Yukawa like potentials with different spin and isospin dependence. The 
simple increase of the K coupling does not solve the problem by itself. For instance, 
for Ak '^Ak and Bk "^Bk, the contribution of all three pseudo-scalar mesons is 
(when b = 1.51 fm): F^m = 1-404 r°, Fp = 0.815 r° and F^/p = 0.723, and when the 
vector mesons are added one gets: T^m = 1-714 F", Tp = 1.130 F° and Tn/p = 0.518. 
Namely, Fjvm turns out to be too large. But from the previous discussion, in relation to 
equations (5.4), (5.5) and (5.6), we have learned that it could be possible to reproduce 
at the same time the data for all three observables by: (i) making the total tensor 
interaction in F^^ small, and simultaneously (ii) decreasing F^^ and increasing F^^, 
without modifying Fnm too much. The first goal can be accomplished, for instance, 
through the modifications: 3i3^ and Bk — ^ '2Bk, and the second one with 

Ar, — > SAr) and Aki —>■ 5Aki- The following cases are illustrated in Table V: 

• {PS'): Only the pseudoscalar mesons are included with the above changes in r] 
and K meson couplings. 

• {PS' — K*): The K* meson potential is incorporated but with the inverted sign. 

• {PS' — V) : All vector meson potentials are included with the inverted signs. 

No best fit to data has been attempted. Yet, it is clear that there are many other set of 



parameters that reproduce reasonable well the data. We wish to stress as well that, when 
the vector mesons are considered, the correct values of Tn/p are obtained only by overturning 
the signs of the vector meson potentials. 



VI. SUMMARY AND CONCLUSIONS 



A novel shell model formalism for the nonmesonic weak decay of the hypernuclei has been 
developed. It involves a partial wave expansion of the emitted nucleon waves and preserves 
naturally the antisymmetrization between the escaping particles and the residual core. The 
general expression (2.13) is valid for any nuclear model and it shows that the NM transition 
rates should depend, in principle, on both: (i) the weak transition potential, through the 
elementary transition amplitudes M.{pPlLXSJT; j\j]s[,mtj^), and (ii) the nuclear structure, 
through the two-particle A'"A parentage coefficients {Ji\ \ (ojjvmtj^ "^Ja) j 1 1 -/f) ■ The explicit 
evaluation of the matrix elements Ai is illustrated as well. 

Two nuclear models for even-mass hypernuclei, namely the EPHM and the QTDA, were 
worked out in detail, and the Eqs. (4.4) and (4.7) were derived. The last one explicitly 
depends on the initial and final wave functions. But, because of: i) the inclusive nature 
of the nonmesonic decay, and ii) the peculiar properties of the coefficients F^^^j{ja) this 
dependence is totally washed out for all practical purposes. In this way we have arrived at 
a very simple result for transition rates, given by the Eq. (4.9), which except for the BCS 
pairing factors Vj^, agrees with the EPHM result. Thus, it can be stated that the two-particle 
correlations in the initial and final state are only of minor importance if of any. With some 
additional effort can also be incorporated the higher order nuclear structure effects such as 
the four quasi-particle excitations, collective vibrations, rotations, etc.. Yet, it is hard to 
imagine a scenario where the later could be relevant at the same time that the former are 
not. Therefore, we conclude that the nuclear structure manifests basically through the factor 
^mt^j{ja)i which is engendered by the Pauli principle; ja stands for the hyperon partner in 
the initial state, and runs over all proton and neutron occupied states in the initial nucleus. 
It is amazing to notice that the Eq. (4.9) is valid for any even-mass system, which can be 
so light as \H and \He are or so heavy as ^°®P6 is. (A quite similar result is also obtained 
for the odd-mass hypernuclei, and this issue will be discussed elsewhere.) One should also 
add that the last equation contains the same physics as the Eq. (5) in Ref. [13] or the Eq. 



(30) in Ref . [27] , with the advantage that we do not have to deal with spectroscopic factors. 
Of course, neither the initial and final wave functions are needed. 

Attention has been given to the nonrelativistic approximation, used to derive the weak 
effective hypernuclear one-meson exchange potentials (3.14) and (3.15). Same errors and 
misprints that appear in the recent papers [13,19,24] have been corrected. Additional parity 
violating vector meson operators <Tjv • f^^\r) and a a ■ i^~^\r), usually neglected, have been 
considered as well. The matrix elements of these new terms were fully discussed, and it was 
found that they are quite important quantitatively and therefore should not be omitted. 

With the OMEM parametrization from the literature [13], and keeping the treatment 
of the FSI at the simple Jastrow like level {g{r) — 1 — jo{qcf)): we reproduce satisfactorily 
the data for the total transition rate {V^nm — r°), but the n/p-ratio {Tl^/p ~ 0.42) and the 
proton partial width (Tf > O.GOPO) are not well accounted for. More elaborate treatments 
of the FSI increase sensibly the n/p-ratio, but they are unable to solve the puzzle [13,24], 
especially after the last experimental result for this observable [43] . We have found that the 
new vector meson operators are not of much helpful in this regard either. 

Finally, bearing in mind the phenomenological nature of the OMEM, we have also tried 
to reproduce all three data simultaneously by varying the coupling strengths in a significant 
way. As the only guide were used the simple formulas (5.6), (5.7) and (5.8), which come out 
from the multipole expansion done within the EPHM. Such an attempt was successful, and 
we get: 0.87 <r*^^/r°< 1.30, 1.06<r^')p< 1.42, and 0.41 <rtVr°^ 0.55. We are conscious 
that changing a coupling by up to a factor of 5, with the sole justification of accounting for 
the data, is a rather a desperate way out of the Tn/p puzzle. Thought no profound physical 
significance is attached to the "new" parameters, it even can be said that such a procedure 
is not physical. However, after having acquired full control of the nuclear structure involved 
in the process, and after having convinced ourselves that the nuclear structure correlations 
can not play a crucial role, we firmly believe that the currently used OMEM should be 
radically changed. Either its parametrization has to be modified or additional degrees of 
freedom have to be incorporated, such as the correlated 27r from Ref. [27] or the 4-barion 



point interaction from Ref. [26], avoiding clearly the double counting. In fact, it would be 
very nice to see the outcomes of such studies. 
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APPENDIX: NUCLEAR MATRIX ELEMENTS 

Here we evaluate the transition matrix elements that appear in (2.15) for the potentials 
V(r, S1S2) defined in (3.14) and (3.15). 

The PC potential contains the operators fuif), /m(^)(*''a ' ^n), and fM{r)SAN{r), and 
the corresponding matrix elements read: 

{pPlLXSJ\fM{r)\mLX'S'J) = SuSxySss'SLL{Pmi){pl\fM{r)\l\), 
{pPlLXSJ\fi,{r){(TA ■ (rN)\miX'S'J) = Su5xx'Sss'SLL{Pmi){pl\fi,{r)\l\)[2S{S + 1) - 3], 
(pPZLA5J|/^(r)5AAr(r)|llNLA',S'J) = (-)^+'+'^+^555'(^LL55iV^AA'[(PL|NL) 



X {pl\fl,{r)\l\) < 


X' 1 J 




/ 1 2 




> < 


> 




1 X 2 




A' A L 



> (/020|I0), 



(Al) 



with 



(PL|NL) = Sll J R''dRjL{PR)nuL{R), 



and 



(p^|/m|1I) = y r^drji{pr)fM{r)Un{r),etc. 
The PV potentials are of the form 



V^^(r, si, S2) ~ S • fW(r) with S = < 



ICTa X <Tn 



(A2) 



(A3) 



(A4) 



and we obtain: 



{pPlLlSJ\V^^\nmU'S'J') = SlJI'\{\010\10) 



1 1 L r 




r s' j] 




II 











X (-)^+^+'+^(5||SM||5')(P^I/l?^|nl)(Pi^|NL). 

The spin dependent matrix elements are: 

1/2 1/2 S' 



(SIMS') ^Vess'i-Y 



1 S 1/2 



i-f^^' (SIMS'), 



and 



(S\\t{cTj, X o-;v)||^') = Vl2(5so55'i + <55i<55'o) = -{SMcTr, X cta)||^')- 
The matrix elements (p^|/M^|nl) are easily evaluated and one obtains, 

(pZ|/l^)|nl) = (P^I/M|nl)±(pi|/i^|nl). 



with 



0>i|/«|nl) ^ / r'dr3,(pr)fM{r) (i^ 



and 



Note that the "sum rule" 



(p/|/M|nl) - {M\n\) = {pl\f'M\nl) = J r'drji{pr)f'M{r)n.^{r), 

should always be obeyed. 

The radial integral (AlO) can be expressed as 



(P'l/M|nl) = -P j r'dr1l„,{r)fM(r) 



2(2Z + 1) 



2(2Z + 1) 



which immediately leads to 

(p^I/mIhI) = Tp/ r^drMpr)fMir)n,i{ry, for Z = I ± 1. (A13) 
We are interested here only in 



and, in order to simplify the integral (A9), the following relationship can be used 
We obtain 

J r^drji{pr)fM{r)Tlu{r); for / = I + 1 , 



(A16) 



/(r^ — \'^b'^)rdrji{pr) fM{r)TZi\{r); for Z = I — 1 . 

It should be remembered that the radial wave functions TZn\{r) and TZul{R) have to 
be evaluated with harmonic oscillator parameters \/2b and b/ \/2, respectively, being b the 
oscillator length for the harmonic mean field potential. 

As indicated in (5.1) the FNS effects are incorporated directly in the radial integrals 
through the replacements /m(^") — ^ /m(^)) ^^c. At variance, the SRC, given by (5.2), are 
added by the substitutions 



|llm) — »■ |llm) = gi{r)\l\m)] \plm) — >■ \plm) = gf{r)\plm), (A17) 

in (Al) and (A3), and when the FNS and the SRC are included simultaneously, the radial 
integrals (A3) become: 

(p^I/m|1I) = / r'drji{pr)gf{r)fM{r)g,{r)n,,{r). (A18) 

Thus, it is equivalent to comprise the SRC either through the wave functions, as done in 
(A17), or by renormalizing the radial form factor: /m(^) ~^ 9f{''^)fM{''^)9i{''^)- The same is 



valid for /^(r), /^(r) and fliir). On the contrary, for the integrals (A9) and (AlO) which 
contain derivatives, from (A17) one has: 



(p^I/mIII) = - r dr\ji{pr)g'f{r) ± pj\{pr)gf{r)]f M{r)gi{r)nu{r); for / = I ± 1, 



Ir'drji{pr)gf{r)f^{rM(r) - rb-Vi(r)]7ei,(r); for Z = I + 1 , 



J r^drji{pr)gf{r)f M{r)[g'i{r) - (rb ^ - |V ^)gi{r)]nn{r); for / = I - 1 , 



(A19) 

being g'{r) = dg{r)/dr. In this case it is no longer possible to include the SRC via the 
form factor, which is a direct consequence of the fundamental difference between the FNS 
effects and the SRC. Namely, while the SRC modify the nuclear wave functions, the FNS 
renormalization is done directly on the vertices of the Feynman diagrams that determine 
the one-meson exchange transition potential. 

Finally, the isospin matrix elements needed in the calculation are: 

(1, -lITol-i -\) = 1; (1, -l|Ti|-i, -I) = 1, 



(0,0|7^|-i, I) — (1,0|7^|-|, I) — 



(A20) 
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TABLE I. Isoscalar (r = 0) and isovector (r = 1) coupling constants in units of 
Gprnl = 2.21 x 10"^. 
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Am A'm B'm 
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K* 
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TABLE II. Geometrical factors jIF^^ ^{ja) 
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TABLE III. Parameters used in the calculations: masses (in MeV), cut-ofFs (in GeV) and the 

isoscalar (r = 0) and isovector (r = 1) coupling constants (in units of 10^^^ MeV~^). 
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T<' 
J\ 
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1.2 


-0.828 




0.228 
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K* 
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1.016 




r = 1 


vr 
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-0.546 




K 
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1.2 


-0.127 




0.764 




P 
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1.4 


0.273 


0.105 


-0.907 


-0.407 


K* 


892.4 


2.2 


0.514 


0.324 


1.072 


0.274 



TABLE IV. Parity conserving (PC) and parity violating (PV) nonmesonic decay rates for j^^C, 
in units of = 2.50 • 10~^ eV. All coupling constants and the cutoff parameters are from Table 
III and b = 1.51 fm. All calculations were done within the EPHM, except for a few results which 
were evaluated in the QTDA and are shown parenthetically. 
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(0.250) 
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TABLE V. Parity conserving (PC) and parity violating (PV) nonmesonic decay rates for ^^C, 
in units of r° = 2.50 • 10~^ eV. The data are taken from Refs. [39-43], and large experimental 
errors are due to the low efficiencies and large backgrounds in neutron detection. The calculations 
were performed for both h = 1.51 fm and h = 1.75 fm, being the later given parenthetically. 
In Calculation A all parameters are from Table III, and PS and V stand, respectively, for the 
pseudo-scalar (vr + 77 + K) and the vector [p + lo + K*) mesons, while the label (P) indicates that 
only the moments Pm are considered (see (4.17)). In Calculation B the coupling constants listed 
in Table III are modified as: ^ 3^^, Ak^ 5^^^, SBj, and Bk ^Bk, and the signs 

of all vector meson potentials are inverted. 
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Calculation A: 
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Calculation B: 








PS' 
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